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Abstract 

Algebraic solution of RPA equations for nucleon re-interactions in 
the case of quasi-elastic charged current neutrino-nucleus scattering 
is presented. Abelian algebra of matrices allows to extract four inde- 
pendent corrections to cross section separately. Results of numerical 
computations are shown. 

1 Introduction 

A better theoretical understanding of nuclear effects in neutrino-nucleus 
scattering is important in view of data analysis from new more precise 
neutrino experiments such as K2K, MINOS, MiniBoone. Determi- 
nation of parameters governing neutrino oscillation phenomenon, in 
particular of ^13 requires an improved knowledge of neutrino-nucleus 
cross sections In the above mentioned experiments neutrino beam 
energy is of the order of 1 GeV. A characteristic feature of neutrino- 
nucleus reactions in this energy domain is formation of resonances and 
subsequent pion production. However, the quasi-elastic contribution 
is still important and its precise determination is of interest. 

Nuclear effects in neutrino-nucleus interaction are often evaluated 
in a framework of Monte Carlo approach. Scattering is split into two 
steps. Neutrino interacts first with a free nucleon, outgoing particles 
are then subject to re-interactions inside nucleus. In more system- 
atic theoretical approaches Mean Field Theory with relativistic Fermi 
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gas of protons and neutrons as a ground state can serve as one of 
techniques to provide a model for nucleus [2|. In its simplest version 
nucleus forms a sea of fermions with momentum uniformly distributed 
inside the Fermi sphere. The effect of MFT is that one has to substi- 
tute nucleon mass M by an effective mass M*. It is well known that 
for energies in the GeV range in the case of electron-nucleus scatter- 
ing Fermi gas model with fine tuned values of Fermi momentum and 
effective mass accounts for basic features of the dynamics 0. More 
realistically, nucleons interact with each other exchanging pions and 
p mesons and also short range correlations have to be considered by 
introducing suitable contact interactions terms In the ring approx- 
imation ot the RPA approach a summation over all Feynman diagrams 
is substituted by a sum of diagrams, where only Ip-lh (one particle - 
one hole) excitation are included In order to make the theory bet- 
ter one should also consider elementary 2p-2h excitations in order to 
enlarge cross section in the so-called "dip" region [fi,J[j. It is however 
difficult to include this contribution in the RPA scheme jH]. 

In this paper analytic expressions for four contributions to RPA 
corrections are derived in the case of quasi-elastic neutrino reactions. 
An algebra of matrices is introduced to solve Dyson equation. Results 
for separate contributions are presented. It is known jH] that RPA 
corrections typically reduce the maximum in the energy transfer dif- 
ferential cross section by a factor of about 10%. Our formalism when 
applied to the CC process is only expected to reproduce these results. 
In particular we want to mention here the paper ^H]- We try to keep 
the same notation in order to make comparison easier. A first motiva- 
tion for the present study is to construct a general framework in which 
more detailed analysis of quasi-elastic CC processes could be possible. 
One can evaluate significance of uncertainties in various parameters: 
nucleons form factors, coupling constants, effective mass etc. The sec- 
ond motivation is that the same algebraic framework can be applied to 
NC reactions and hopefully also (with necessary modifications) A ex- 
citation. In our RPA computations we keep a constant value of Fermi 
momentum. Also M* is assumed to be a function of Fermi momen- 
tum only thus it is a constant. Inclusion of local density effects in the 
analytical framework is in principle possible but rather complicated - 
a lot of manipulations with spherical harmonics are necessary [Jj with 
approximations difficult to control. For all practical purposes it is suf- 
ficient to have exact cross section formula for a fixed value of Fermi 
momentum (and M*) since one can perform a numerical integration 
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over Fermi momenta with a distribution defined by a density profile 
of a nucleus in question. Our algebraic scheme enable rather simple 
computation of both effects (RPA and local density) and this is the 
third motivation for our work. We wish to mention that similar idea 
of solving Dyson equation can be found in 

The paper is organized as follows. A short description of the model, 
Feynman rules and the main technical tool - algebra of matrices is given 
in section 2. Explicit expressions for RPA corrections are presented in 
section 3. Section 4 contains a discussion of results and also a compar- 
ison with RPA computations done in the relativistic generalization of 
the Marteau model |T^. There is a very nice agreement between two 
approaches. Some technical details of algebraic computations are col- 
lected in Appendix. The aim of this paper is to present main features 
of the approach. More details and discussion will be contained in the 
next paper being in preparation. 

2 The formalism 

We consider charged-current (CC) quasi-elastic neutrino-nucleus scat- 
tering (Fig. P). The model of nucleus is given by the Mean Field 

[i{k') p(p') ^l (k') p(p') 




v(k) n(p) v(k) n(p) 

Figure 1: The basic diagram describing neutrino-nucleus interaction. Scattering takes 
place on a single nucleon with a definite momentum given by a Fermi gas distribution. In 
energy domain of few GeV effective four-Fermion vertex provides an excellent approxima- 
tion. 

Theory j2| (nucleon mass becomes an effective mass M*{kF)) with in- 
teractions due to contact terms and exchange of pions and p mesons 
0. In a first approximation the nucleus is treated as a relativistic 
Fermi gas with Fermi momentum kp determined by nucleons density. 
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Elementary weak charge current nucleon-nucleon current is ex- 
pressed by means of form factors 



We omit Gp{q^) term since its contribution to and z/^ cross sections 
at El, ~ IGeV is negligible. It can be put into our scheme if required 
with only minor modifications. 

The differential cross section (per nucleon) reads: 



d\q\ dqo IQn'^ppE'^ 

Pf = kp/Sn'^, we assume the same values of Fermi momentum for 
neutrons and protons. L'^'^ is a leptonic tensor: 

L^u — 8 (jt^k I, ^ k ^ky — g^yk^k i 'i^^uapk k^^ . (3) 

The sign ± depends on a process considered (neutrino/antineutrino). 

The polarization tensor H^^ is a basic object that contains full in- 
formation about nuclear effects. It is defined as a chronological product 
of many body currents: 

W^'igo, q) = -^J d'xe^'?""" < 0|T {J^'{x)J''{0)) |0 > . (4) 

The |0 > is a ground state of nucleus described by the Fermi gas 
model. Prescription for nuclear physics [Tij enables evaluation of U'^'^ 
by means of the standard QFT techniques with modified (depending 
on kp) progagator G{p). 



2.1 Feynman rules 

The polarization tensor is split into a "free" part and RPA correction: 

U>^^ = U>}:^^ + AU^^^p^. (5) 

The "free" tensor is given by a simple fermion loop (Fig. |2l which is 
spanned between two vertices with form factors insertions. 

^Treeil) = J T^lTr {G{p + q)r^ {q)G{p)r'' {-q)) (6) 



4 



r" 



Figure 2: li^^^^ 

G{p) = {k+ M*) + ^<5(po - E,)eikF - p)j . (7) 

G{p) describes a propagation of a free fermion in the Fermi sea. 

RPA corrections arise from summation of an infinite sum of Ip-lh 
contributions p]. These corrections are given by Feynman diagrams 
containing tt propagator - V^'^ , p propagator - W^'^ and interactions 
vertices NNn as well as NNp. Landau-Migdal parameter g' is put 
together with genuine pion propagator to form a redefined pion "prop- 
agator" (for details see [El). 

In what follows in this section we do not write down explicitly 
Lorentz indices. It will be understood that unless specified all the 
objects are 4x4 matrices with indices ( )^'^ . We define tensors Hp 
and IItt as loop diagrams with CC weak nucleon-nucleon current in one 
vertex and NNp or NNtt vertices (FigOI). Notice that contributions 
from tensors given by figures (jHk) (|Ht>) are equal. The same applies 
to contributions ([Sh) and ^^). This property simplifies the algebraic 
form of RPA corrections. 

'^RPA is an 8 X 8 matrix defined by an infinite series: 

Arpa = Ao + AoHgAo + AonGAoncAo + (8) 

which is illustrated in Fig. |^ We defined two new tensors -8x8 
matrices: Hq and Aq. Their definitions can be simply understood in 
diagrammatical language (Fig. 

With all these definitions the formula ^ can be rewritten in the 
form of Dyson equation (Fig. d): 

Arpa = Ao + AqUgArpa- (9) 
It is clear that RPA corrections to polarization propagator are given 
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a; 



n = 




( K propagator subtracted ) 




c) 



n = 



d) 



( p propagator subtracted ) 





Figure 3: Diagrammatical definitions of Uj^^p. They are loops with CC weak nucleon- 
nucleon current in one vertex and NNtt or NNp vertices. 



-<m- 
m>- 



-<m 



Figure 4: Following |^ we define RPA propagator as a sum over Ip-lh diagrams with 
external meson fields propagators. 



Figure 5: The propagator RPA can be expressed as an infinite sum. 
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Figure 7: Dyson equation for RPA propagator. 



AIIrpa = 




Figure 8: Diagrammatical explanation of RPA corrections to polarization propagator. 
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by ^RPA multiplied from both sides by Hp and 11. 




The diagrammatic explanation of this formula is presented in Fig. 
IHl Our strategy is to solve equation for /^rpa and then to obtain 
expression for AUrpa- 



2.2 Algebraic properties of the polarization 
tensor 

A coordinate system is chosen in which four momentum transfer reads: 



9m 



kp - k'f, = (go, 9,0,0). 



fill 



The Dyson equation is a 8 x 8 matrix equation. In order to solve it we 



introduce 4 x 


4 matrices: 








/ 




qoq 





o\ 




yfi 








CL = 




4 


























V 










0/ 






/I 








'\ 









1 








CA = 








1 









\o 








1/ 



ex 



/O 0\ 



10 

Vo i/ 



evA 



They satisfy matrix multiplication relations: 



/O 





Vo i 





—i 
/ 



(12) 





= ex 




ei CT 


= = 


ex ex 


CL eA 


= ex = 


eA ex 


e-L evA 


= = 


evA ex 


er ex 


= ex 




ex eA 


= ex = 


eA ex 


ex evA 


= evA = 


evA ex 


GA eA 


= eA 




GA evA 


= evA = 


evA eA 


evA evA 


= ex- 





(13) 
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The polarization tensor must be a linear combinations of ex, ct, 

n = u^eL + n^cT + n^^evA + n^e^. (14) 

This is clear since IT is given by a sum over Feynman diagrams 
and each contribution is of this form due to the fact that algebra 
eL,eT, evA, &A is closed under multiplication and all the building blocks 
are expressed in terms of these four basic matrices. Consequently the 
cross section has the form 

(15) 

where Ll = Ly^ei^^ etc. Many authors, e.g. ^01 call four contribu- 
tions: longitudinal, transverse, V-A, and axial. This can cause some 
confusion because contributions to the cross section are sometimes 
called after spin-isospin operators present in nucleon-nucleon transi- 
tion current. We will return to this point in Section 4. 



3 RPA corrections 

A lot of simplifications come from the fact that Aq contains only lon- 
gitudinal, transverse and axial terms |lfl| . 

. _(W \ _ fW^eL + W^CT + W^CA 

(16) ■ 

W and V are 



W^= = — o — (17) 

qf,-mf, + le 

. QT, — m,, 1 

= „ % V— (18) 

— -|- le 

= = - , . (19) 

= , . -g' (20) 



(21) 



g' is the Landau Migdal parameter. 
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Also 

Up = n^ei + nJeT + n^^evA (22) 
= U^eL + Uler + uVevA + ^ieA (23) 

We find 

- T-T T-r - rtVA„__ . ttL „ . , ttT „_ , ttA „ . \^^) 



We introduce a general notation: 

= Af CL + Af eT + ^VevA + ^Ua i = 1, 2, 3, 4. 

The 8x8 matrix equation ^ can be rewritten as a set of four 4x4 
matrix equations: 

Ai = + T^n^pAi + wn^^Ag 
As = W^HppAa + T^np^A4 

A3 = yn^^Ai + yn^^A3 

A4 = y + yn^^As + yn^^A4. (26) 

Equations (|26)l are transformed to a set of algebraic equations and 
solved in the Appendix. We obtain a general solution for ^rpa 

_ /A^^e^ + AfeL + AfcT A^^eyA \ 

^^^^"l ^V^YA Afe^ + Afei + Alerj- ^^^^ 

RPA corrections to polarization tensor are: 

An^p^ = ((n^)2 + 2n^n^) (Af + Af) 

+ (n;^)2Ai + (n^)2 (a^i + Af ) (28) 
An^PA = ((n^)2 + 2n^n^;^ + (nr^)2) (Af + Al) + (n^)2Al 



(( 



nr)^ + «-^)2 Af + Ai^ 



+ (n^^-^n^ + u^^ui + njnr ) (a^ + Ar^^) (29) 

An^#A = (nJn^ + nX + nrnr)K^ + A3^^) (30) 

+2 (Af + njf^nj + 2 (aI + Af ) (n^^ + n^) n^^ (31) 

An^P^ = {U^f^i (32) 
Expressions for A^... are written in the appendix. 
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4 Discussion 



Numerical evaluation of RPA corrections to H^'^ requires a knowledge 
of both real and imaginary parts of all objects presented in (|28l - [82)l . 
We get the necessary formulae from ^H]- With this input one can 
perform numerical analysis of our algebraic results. 

We present some plots demonstrating consistence of our proce- 
dures. We compare RPA corrections obtained in our paper with those 
computed in a different approach. The approach we choose is based 
on Marteau model with some improvements A use of relativistic 
generalization of the Lindhard function made kinematical regions of 
both models the same. 

In the Marteau model three contributions to cross section are iden- 
tified according to spin-isospin operators present in the transition am- 
plitude. To find a bridge between two decompositions we observe 
that hadronic tensor components H^^ , H^^ , H^^ , H^^ contribute only 
to longitudinal and charge contributions (in spin-isospin nomencla- 
ture) while remaining components contribute only to the transverse 
part. This is strictly speaking true in the approximation when (pis 
target nucleon momentum) terms are neglected in H^^'^ which is valid 
within few %. We decided therefore to single out two contributions 
in both approaches and to call them in order to avoid confusion: I 
and II. The contribution II is equivalent to the sum of charge and 
longitudinal parts in the Marteau approach while the contribution I is 
equivalent to the transverse part. 

The identification of I and II parts in our approach requires some 
algebra. We obtain: 

(Vn^i =Lt{u^ -U^)+LvaII''^ (33) 
(L^.n^'^)!! = LlU^ + {La + Lt) (34) 

In numerical calculations we assumed the following values of param- 
eters present in the theory: Ma = 1.03 GeV - axial mass, standard 
values of coupling constants for pions and p mesons, Landau Migdal 
parameter g' = 0.7 except from two comparison plots (in jT2] the value 
g' = 0.6 was assumed and we take the same value in order to make 
the comparison consistent). Effective mass was calculated according to 
self-consistency equation of MFT theory [2]. It is assumed that target 
nucleus is oxygen ^^O and that Fermi momentum is kp = 225 MeV. 
We get M* = 638 MeV. 
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CC quasielastic neutrino total cross section on sO''^ 




log(E„[GeVl) 

Figure 9: Comparison of the computations of quasielastic neutrino total cross sections 
on oxygen, g' = 0.7. 




Figure 10: Comparison of differential quasielastic neutrino cross section for Fermi Gas 
of free nucleons with cross section modified by RPA corrections for neutrino energy 1 GeV. 
The calculation was done for the free nucleon mass and g' = 0.7. 
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CC quasielastic differential neutrino cross section on gO^® with IVI = 638 IVIeV 




900 



Transfer of energy: qg [MeV] 

Figure 1 1 : Comparison of differential quasielastic neutrino cross section for Fermi Gas 
of free nucleons with cross section modified by RPA corrections for neutrino energy 1 GeV. 
Calculations were done with effective mass M* = 638 MeV, g' = 0.7. 



CC quasieiastic neutrino differentiai cross section on gO^^ witfi U = 939 MeV 




too 200 300 400 500 600 700 800 900 

Transfer of Energy qg [MeV] 



Figure 12: Comparison of differential quasielastic neutrino cross sections with RPA 
corrections in: the model of this paper and Marteau modelT3| for neutrino energy 1 GeV. 
The calculation was done for the free nucleon mass and g' — 0.6. 



13 



CC quasielastic differential cross sections on gO with M = 939 MeV 



Part 11, Marteau Model 
Part I, Marteau Model 
Part II, Computations of this paper 
Part I, Computations of this paper 




300 400 500 600 

Transfer of Energy qg [MeV] 



Figure 13: Comparison of RPA corrections in contributions I and II, in the model of 
this paper and Marteau model U2|. The calculation was done for the free nucleon mass, 
- 0.6. 



In the FigEl we compare predictions for the total cross section in 
three cases: (i) free Fermi gas with M* = 939 MeV, (ii) RPA compu- 
tations with M* = 939 MeV, (iii) RPA computations with M* = 638 
MeV. Inclusion of RPA correlations makes the cross section smaller. 
In the third case reduction of the cross section is more significant for 
neutrino energies up to about 3 GeV. 

In the Fig^land Fig^Jwe show differential cross sections in en- 
ergy transfer for neutrino energy of 1 GeV. As above we distinguish 
two cases in which effective mass is taken either as a free mass of nu- 
cleon or as 638 MeV. One can see the typical expected behavior: in 
the RPA case quasi-elastic peak becomes reduced but at larger values 
of energy transfer an effect of RPA is to increase slightly the cross 
section. We notice that due to effective mass kinematically allowed 
regions in energy transfer are in two cases different. 

In last two figures we compare our differential cross sections with 
predictions of the model described in ^J]- A good agreement between 
influence of RPA corrections in two models is seen. The contribution 
I is dominant in both cases. Differences between them are small. The 
behavior of contribution II in both cases is similar. Marteau model 
gives rise to smaller contributions at energy transfer of ~ 50 MeV and 
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the whole contribution becomes reduced by about 25%. 

We conclude that our algebraic solution of RPA equations leads 
to modifications of cross section similar to other approaches. We can- 
not expect that Marteau model can produce numerically identical 
results as it is a hybrid model which combines nonrelativistic poten- 
tial approach with relativistic Lindhard function. We hope that our 
algebraic scheme will be useful in other cases mentioned in the Intro- 
duction. 

5 Appendix A 

Equations are rewritten as sixteen: 

At = 

= 

At = V^Aiut 

At = y-^(i + A^n^,) (35) 

Af = (^Af + Af ) (W^ + ly^) U^p 

A^ = (a^ + A^) (w^^ + VF^) n^^ 
Af = AiiutV^ + UtV^+UtV^) 

+Af {v^ + v^) (ni + n^,) 
Ai = y^Af (nlyVntF^ + n^y^) 

+A^{v^ + V'^)(ut+Ilt) (36) 

Af = + (a^^ + Af ) nj^ (vF^ + py^) 

+A3^X-^(^^ + ^'^) 
Ai^ = (A^ + A^)nJ,(w^ + Ty^)+ArX/(^^ + ^^) 

= ArXV(y^ + y^)+Al^(nl + nL)(i/-^ + F^) 
+At{utv^ + ul^v^ + ul^v^) 
Al = y^ + ArX/(y^ + y^)+Al(nL + nl)(y^ + y^) 
+At (utv^ + ul^v^ + n^v^^) (37) 



15 



{A^nl + Af n^,^ + Aln^V) {w^ + 
{Atu^^^ + Af n^,^ + A^^nl + A3^^nL) {v^ + 1^^) 
((A^ + A^n^,^ + Ar^ni + A^nL) {v^ + v^) (38) 

We solve these equations sector after sector. Equations ([SSI -ESI lead 
clearly to: 

. A _ tjtA 



^4 



A^ 
Af 








1 _ yAjjA 

To proceed further it is convenient to define: 



(39) 



R 

V' 
V 



TA 



TA 



LA 



(40) 

Equations (|36l - l36|) contain only A:^^ and A^ components. We obtain: 



Af 



A^ 

^4 



1 - RLAULp 








yL I yLA-^L yi 



(1 - v^nt) (1 - v^^{u^^ + ut)) ■ ^^^^ 

A'^ and A^^ components mix among themselves but always in pairs: 



A^ 

Af 
Al 



/\VA 
^3 

Ar^ 

aP 
aX^. 



(42) 



We derive: 



Af 





>-^^^(nl + nL) 








[i-y^^(n4 + n^J] 




- i?^^y'r^(n)f/)2 
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^rAJ^TAyTA^-^VAy 



[1-VTA(UA^ + UT^)] 



_ JlTAvTA^YlVA)2 



(43) 



[1 - VTA (nA^ + n^J] [l _ i^TAnT ] _ yTA^TA(nVA) 



^AYTAjfAmVA\2 



[1 - (n^, + n^j] 



RTAUT^ 



- V^AjiTA^iiVAy 



(44) 



A^^ 



^3 



A^ 



= 



w^^n^/(y^ + (i-n,,y^)Af^ 



1 - W^^(nj, + VAT(UVAy^ + VAT(UA^ + UT^){UT^WAT if 



1 - w^^(nj, + F^^(n^^)2) + FAT7nA^ + ut^){ut^wat l^Q 
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